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At the BCS level of chiral symmetry breaking, the mass 
gap equation for quark-antiquark condensation only uses the 
kernel of the Bethe-Salpeter equation. We introduce coupled 
channels with ladder mesons in the mass gap equation. Con- 
sistency is insured by the Ward Identity for axial currents, 
and the 7r remains a Goldstone boson in the chiral limit. We 
find that bare mesons with confined quarks do not contribute 
to the mass gap equation, and estimate that the contribution 
of full mesons may be large. 
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The studying of dynamical symmetry breaking beyond 
BCS has never been done before. This amounts to join 
the BCS mechanism with the mean field expansion of ef- 
fective bound states. Whenever a field theory produces 
bound state solutions, one should try to include them at 
the onset in the self-consistent Schwinger Dyson equa- 
tions of that theory. The solution of this technical prob- 
lem is not only relevant from condensed matter physics 
to particle physics, but it is essential for a correct under- 
standing of hadronic physics. 

Strong interactions were observed, at the very begin- 
ning of hadronic physics, with large scattering lengths 
and wide decay widths. The discovery of quarks as the 
microscopic elements of hadrons, and two new theoret- 
ical concepts, confinement and chiral symmetry break- 
ing, eclipsed partially the importance of coupled chan- 
nels in theoretical hadronic physics. Several different 
strong or confining, chiral invariant quark-antiquark in- 
teractions have been used to condense in the vacuum 
scalar quark-antiquark pairs, inducing dynamical chiral 
symmetry breaking, and a Goldstone boson, the tt was 
found Q, in the case of vanishing quark masses. Re- 
cently coupled channel effects have been reevaluated in 
many hadronic phenomena. A paradigmatic case is the 
omnipresent nucleon, where the estimated negative mass 
shift due to the coupling to the channels tt N, tt A, 
tt N * . . . has increased almost an order of magnitude 
in the last decade Moreover in the literature only 

1 or 2 coupled channels are usually included, and it re- 
mains to be proved theoretically that the total mass shift 
due to the infinite series of coupled channels is finite. 

The crucial problem seems to dwell in the pion mass 
. The pion itself is coupled to the channels tt p, p p- ■ ■ 
and his bare mass is believed to suffer a large shift of the 
GeV order. This effect is so large (|] , that it disputes with 
the instanton S the solution of the U(l) problem of the 



tt rj mass shift. Now suppose that the mass gap equation 
for chiral symmetry breaking was solved at the BCS level, 
i.e. without including the coupled channels, then a bare 
pion with vanishing bare mass would be a solution of the 
Bethe Salpeter equation. If the coupled channels were 
then included, at posteriori, in the bound state equation 
then the pion mass would be the sum of the small bare 
mass plus a large negative mass shift and thus would have 
a resulting negative mass, which is absurd. We will have 
to use the Ward identities (WI) in order to insure that 
the bound state equation for the pion -a Bethe Salpeter 
equation with coupled channels - is consistent with the 
non linear mass gap equation. A reward of this study is a 
tt with a vanishing positive mass in the chiral limit, and 
a tower of states above the tt with higher masses [Q . 

The WI were first derived for QED, and concerned the 
free vector vertex T'j — 7 M , and the free Dirac fermion 
propagator <S/(fc) = i/ /k — m + ie. There is also a WI 
identity for the free axial vector vertex I^ 5 = 7^7 5 that 
involves the free pseudoscalar vertex = 7 s , 

-i{k M - /fcp5(fc)r^ 5 (fc - k')S(k') 

+2imS{k)T\k - k')S(k') = 5(fc) 7 5 + j 5 S(k'), (1) 

where the difference in the right member of the equa- 
tion extends the identity to renormalized propagators 
and vertices. The WI are crucial to ensure the consis- 
tency between the mass gap equation (MCE) and Bethe- 
Salpeter (BS) equation. They ensure that the self energy 
of the MGE is obtained from the BS kernel by closing 
the fermion line where the vertex is inserted, and that no 
double counting occurs. Inversely they also ensure that 
the BS kernel is obtained if one inserts the vertex in all 
possible propagators of the self energy. 

This will now be illustrated at the BCS level. In this 
case the mass gap equation and the BS equation are, 
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(2) 
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where the full propagator is denoted by — < — and the 
vertex will be denoted by •=. The effective quark-quark 
2-body interaction, a chiral invariant and color dependent 

interaction, is represented with a dotted line With 

the WI, it is possible to check directly that these eqs. are 
equivalent, both for dressed and for free propagators and 
vertices. Nevertheless we will now illustrate explicitly 
how the free WI recovers the MGE from the BS equation. 
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Up to second order in the potential insertions, the vertex 
obtained from eq. (j^) is, 



t'= = r o+ fc'H[] r o + fc' I ° , ° "~^ r o 



+ fc' [H3 ' r o+ fe ' uS r «+ • • • ( 4 ) 



Substituting the WI in eq([|), we get, 

75 S- 1 (fc') +S-\k) l5 = 75 5 - 1 (fc') +5 - 1 (fc) 75 + 

k q 




+ <]■' p J' + p J' + p %/ + ■■■ ( 5 ) 

where the Dirac 75 was represented in the diagrams with 
an *, the loop variables were omitted and we only showed 
the k or k' dependence of the propagators. In the case of 
a chiral invariant interaction, the vertices anticommutc 
with 75 , which cancels exactly the spurious diagrams de- 
pending on both k and k' , and we recover the propagator 
of eq. (||) up to second order in the potential insertions. 

A key product of the WI is the proof that a pseu- 
doscalar Goldstone boson exists when current quark 
masses vanish. We will follow a variant of Pagel's proof 
H which also yields the BS amplitude itself. When dy- 
namical symmetry breaking occurs, the full propagator is 
renormalized and the selfconsistent equation for the self- 
energy £ has a solution with a finite dynamical quark 
mass, 

S-^k) = S^ik) - £(fc) , j£ = A(k)- /kB(k) . (6) 

If we substitute this propagator in the WI, we find the 
solution for the pseudoscalar vertex T 5 with a vanishing 

r 5 (fc = fc') = A(fc) + m 7 5 (7) 

TO 

which diverges for a vanishing quark mass m and shows 
that the pole of a massless pseudoscalar meson appears 
in the axial vector vertex, with a boundstate amplitude 

A y 7 where f v is a norm. Because of the axial anomaly, 
the identity (0) for flavour singlet currents is not verified 
in this form, and it is necessary to extend (trivially) this 
relation to flavoured currents in order to keep track of 
the Goldstone boson. 

The first step to go beyond BCS is to include a ladder 
exchange potential in the self energy £, which ladder in- 
cludes a series of meson exchange interactions. We also 
include an extra effective interaction, which simplifies the 
the energy dependence of E, and limits the meson cou- 
plings in the bound state interactions to vertices of 3 
mesons. The mass gap equation is now, 



where — »Z — I CZH represents a ladder of effective potential 
insertions in a pair of quark and antiquark lines. Us- 
ing the prescription of inserting the vertex in all possible 
propagators we arrive at the corresponding bound state 
equation for the vertex, 

•= = r + uZ>= + 




where the first line corresponds to the BS equation at 
the BCS level. The last line, which includes the meson 
exchange interactions, has a null effect in flavour vectors. 
It is relevant only for flavour singlets, and may include 
the axial anomaly. 

The quantitative results of this paper will be ob- 
tained with a particular chiral invariant strong potential, 
which is an extended version of the Nambu and Jona- 
Lasinio potential Q. In Hadronic Physics the effective 
interaction should be simultaneously confining, in the 
Minkowsky space, local, and Lorentz invariant, but un- 
fortunately no interaction which complies with all these 
constraints has yet been used to study chiral symmetry 
breaking. Our choice is to relax the last constraint, and 
use an instantaneous approximation which simplifies the 
energy dependence of the interaction. This approxima- 
tion also has the advantage to allow a straightforward 
application to condensed matter physics or low energy 
nuclear physics. We use Q a 2-body potential for Dirac 
quarks, (- § I • § ) ( 7o • 7o ) [Kg (x- y ) 2 - U] 5{t x - t y ) where U 
is an arbitrarily large constant and A are the Gell-Mann 
matrices. The infinite U reappears in the self energy, 
and for color singlet channels this cancels the infinitely 
attractive potential. Any colored state will have a mass 
proportional to U and will thus be confined. The choice 
of an harmonic potential is not crucial, a linear or funnel 
potential have been also used, it simply has a suitable 
Fourier transform — (27r) 3 Afc. Excellent results are ob- 
tained in a very broad range of hadronic phenomena, 
with A'o — 250Af eV, excepting and (^ip) which turn 
out j| to be underestimated by factors of 4.5 and 2 3 . 
With an instantaneous interaction, we prefer to substi- 
tute the Dirac fermions in terms of Weyl fermions, with 
propagators, 

S q (w,k) =S,(w,k) = - E l [k) - ^ (10) 
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which will be represented with an arrow pointing for- 
ward in the case of a quark, and an arrow pointing back- 
ward in the case of an antiquark moving forward in the 
time direction. This formalism is convenient to simplify 
the BS equation into the Salpeter equation, in a form 
which is as close as possible to the simpler Schrodinger 
equation. When Weyl propagators are used, the vertices 
of the effective potential are redefined, they include the 
spinors u and v. The Dirac vertex 70 is now replaced by 
v)u , u^v , v^u or v^v when the vertex is respectively 
connected to a quark, a pair creation, a pair annihilation 
or an antiquark. 

Using the Weyl fermions, we expand the ladder 
in meson propagators k , in Salpeter amplitudes 
4> + , <f>~ and in Bethe-Salpeter truncated amplitudes 
X = fV(u(j) + v^ + v(f)~v)). We find that the self en- 
ergy of the quark has a diagonal component which 
contributes to the dynamical mass of the quark, 




(11) 



and the self energy of the antiquark is exactly the same. 
In eq.(pd|) we only included the nonvanishing diagrams, 
of order U or 1. The BCS quark propagators are propor- 
tional to £7 _1 . The interactions without a pair creation 
or annihilation are proportional to the infinite infrared 
constant U, while the remaining interactions are finite. 
It turns out that the diagrams with no creation or anni- 
hilation vertices are of order U while the diagrams with 
one finite interaction are finite. The self energy also has 
antidiagonal components S a which must vanish for the 
sake of vacuum stability. These components would cither 
create or annihilate scalar mesons in the vacuum. The 
diagrams that contribute to the antidiagonal component 
of the self energy in the mass gap equation are now, 



Sa = 




where the lines . 



(12) 



only mean the presence of a spinor 



u, v, u' or v' and don't include the quark or anti-quark 
Weyl propagator. The first pair of diagrams in eqs. (11), 
( |l2| ) are BCS diagrams. We will now concentrate on the 
new diagrams which depend on the Coupled Channels. 
We will for instance consider the third diagram of eq. 
(|IT|). The integration of the quark propagators, which 
arc functions of the relative energies yields, 



dw 



2tt w — Ti + ie E — w — T2 + ie w — Ty + ie 

i i 



E — T\ —T 2 +ieE — T v - T 2 + ie' 



(13) 



and we obtain a product of propagators Go of quark- 
antiquark pairs. The right Go is in fact part of the ladder 
and should be absorbed in the ket <fi~ . We now integrate 
the remaining Go, together with the propagator of the 
ladder boundstate with energy Mj (P) , 



dE 



2tt -E - Mi + ie E - T x - T 2 + ie 



-Mi -Tx-T 2 +ie 



(14) 



which cancels the potential — iV = (Mi — T\ — T 2 )/i that 
comes from the Salpeter equation (^). After removing 
the poles in the propagators with this integration in the 
loop energies, we get the quark energy T — tp /ku — S f / 
where, 



d 3 k } 

7T 3 

d 3 P- 



d 3 P- 



Sj = " '^—T- + " ( / 2^5> * > 



(15) 



where P is the center of mass momentum of the coupled 
meson. The mass gap equation = iv Jkv — S a with, 



d 3 k rrr uu^ — vir 



/d 3 P \ t t 

~2~~a §~ ) v ^] v + 4>~v^x v + u^xutfi^ 1 - (16) 

Only the small component <p~ contributes to the result, 
and this shows that the effect of the infinite tower of cou- 
pled channels is finite. In fact only the pion yields a size- 
able contribution. The contribution of the remaining infi- 
nite series of mesons, which are excited states and are not 
Goldstone bosons, is negligible because their ^aMf 1 
decrease, and because their average momentum increase. 
In eq. (|15|), the infrared divergent kinetic energy, 



T = rl \ 



d 3 P 
(2tt) 3 



4 1 4> ) + finite terms. (17) 



is now larger than in the BCS case, where it was sim- 
ply £7/2. We will now concentrate on the second line 
of eq. (^) , and analyze the terms proportional to the in- 
frared divergent U. The diagrams with E s are clearly 
in that class, and the last diagram also contributes in a 
subtle way. This diagram can be divided in 4 different 
diagrams, which have a different structure in the vertices 
of the interaction and in the quark propagators. In Dirac 
notation, and skipping the arrows in the propagators for 
simplicity, we can sketch, 
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(18) 



and we find that the first diagram vanishes, the second 
is finite and corresponds to the typical coupled channel 
diagram, and the last two ones are of order U. Bringing 
together the dominant diagrams, we find for the flavour 
vector channels, with the Weyl notation, 




and the dominant terms are, 

Ti + T 2 + T 3 + T 4 - Mi - U 



(Ti +T 2 )G r = U[1+ 2- 



T 3 + T 4 



d 3 P 
(2tt) 3 



G r, (20) 



which is correct at order U, where Mi = T\ + T 2 + (-U)- 
Thus we find a finite solution, with an exact cancelation 
of the terms proportional to U, although the infrared 
divergent part of the kinetic energy is no longer U/2. A 
cancelation also happens in flavour singlet channels. 

We now study the finite effects of coupled channels. 
Using the <f>~ BCS amplitude of the n, we can estimate, 



In this case the (f>~aU~3 is suppressed, and the 
Salpeter equation for the ladder reduces essentially to 
a Schrodinger equation for the (f> + component. The cou- 
pled channel contribution to E a is vanishing, of order 
U~ 3 , and thus the mass gap equation is unchanged. The 
remaining contributions from the coupled channels to the 
bound state equation are of order Mj~ , and disappear 
in the same way. As we predicted the mass of the bare 
(ladder) it meson - which is now proportional to \/U and 
infinite - is exactly canceled by an opposite mass shift 
which is induced by the coupled channels, and the phys- 
ical 7r remains a Goldstone meson in the chiral limit. 

We developed an example of a finite renormalizable 
field theory, and there are good perspectives to apply 
this method to solid state physics. The net result of in- 
cluding the (bare) coupled channels in eqs. is that 
no physical deviation from the BCS case is found. In or- 
der to get the finite contribution from coupled channels 
we must either use a finite U (nonconfining models), or 
replace the bare (ladder) meson exchange in the mass gap 
equation (||) by full meson exchange. It is possible to an- 
ticipate that in this model the effect would be very large, 
in view of the result of eq ( pl| ) which would increase the 
quark condensate by a factor of ~ 2.5 3 . In principle /„• 
might also increase, thus improving the physical results. 

I am very grateful to Prof. J. Emilio Ribeiro for many 
long discussions and, in particular, for pointing the rele- 
vance of the 7r mass problem B . 



d 3 P 
(2tt) 3 



2.57V/, N f 



(21) 



which is large, due Q to the small /„-, but finite. However 
if one goes beyond BCS, the ladder mesons pick up a large 
infinite mass according to Eq. ([l7]) , and this integral 
disappears because ^aMf 1 — > 0. If one then tries to 
solve the mass gap equation iteratively, starting from the 
BCS solution, one finds a pair of accumulation points. 
The first is similar to the BCS solution, and the second 
has a coupled channel contribution which renormalizes 
the interaction, 



T = hj , Mi = 0U 
1 + e 

T = —-U , Mi = eU 



V 2 



(l + 2e) 2 M, 



"772 =0 '(22) 



The stable intermediate solution is close to the BCS one, 
it is obtained when e = (c/U) 2 ' 3 is a function of U, 
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